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1 Introduction 

The development of the fractional calculus dates from the XIX century. Mathematicians as Lacroix, 
Abel, Liouville, Riemann and Letnikov attemped to establish a definition of fractional derivative. But 
the definition given by Caputo in 1967, was an open door to the beginning of the physics applications, 
while the previous definitions enabled a great theoretical development. 

The research on the theory of fractional differential equations has begun to develop recently, and 
in the past decades many authors pointed out that derivatives and integrals of non-integer order 
are very useful in describing the properties of various real-world materials such as polymers or some 
types of non-homogeneous solids. The trend indicates that the new fractional order models are more 
suitable than integer order models previously used, since fractional derivatives give us an excellent tool 
for describing properties of memory and heritage of various materials and processes. Works in this 
direction are P El El m [24]. 

This paper deals with the fractional diffusion equation (here in after FDE), obtained from the 
standard diffusion equation by replacing the first order time-derivative by a fractional derivative of 
order a > 0 in the Caputo’s sense: 

oDfu{x,t) = Uxx{x,t), —oo <x < oo, t > 0, 0 < a < 1, 

where the fractional derivative in the Caputo’s sense of arbitrary order a > 0 is given by 

„ „, , f —s f^(t — n—l<a<n 

[ ’{t), a = n. 

where n G N and T is the Gamma function defined by r(x) = e~^dw. 

Whereas that the one-dimensional heat equation has become the paradigm for the all-embracing 
study of parabolic partial differential equations, linear or nonlinear (see Cannon 0 ), the FDE plays a 
similar role in fractional parabolic operators. 



The FDE has been treated by a number of authors (see [H [H ttH ISl EU) and, among the several 
applications that have been studied, Mainardi |20| studied the application to the theory of linear vis¬ 
coelasticity. 

Generalizations of the maximum principle for initial-boundary-value problems associated to the 
time-fractional diffusion equations were given by Luchko in m and |16| . and uniqueness results there 
were obtained. 

Eberhard Erederich Ferdinand Hopf was an Austrian mathematician who made significant contri¬ 
butions in differential equations, topology and ergodic theory. One of his most important works are 
related to the strong maximum principle for partial differential equations of elliptic type m)- 

In his work |12) . an important theorem related to the sign of the outside directional derivative of 
a function that is a solution to an elliptic partial differential inequality is proved. This theorem was 
proved later for partial differential operators of parabolic type by A. Friedman [7] and R. Viborni |25) 
separately. A weak adaptation of this theorem can be founded in |4], named Hopf’s Lemma, and my 
propose is to generalize it for the FDE. 

That is, under certain conditions that will be enunciated later, if tt is a solution of the following 
problem 


(^) 

oDfu{x,t) = X^Uxx{x,t) 

Slit) < X < S2it), 0<t<T, 0<a<l 

{ii) 

II 

0 <t <T 

(Hi) 

u{s 2 {t),t) = h{t) 

0 <t<T 

(iv) 

u{x,0) = fix) 

a < X < b 


where si and S 2 are given, then, if u assumes its maximum in a boundary point, let us say {s 2 {to),tQ), 
it results that Ux{s 2 {to),to) > 0. 

2 A Fractional Hopf’s Lemma 

Let us consider the moving-boundary problem for the FDE defined in ([1]) where: 

(HI) The curve si is given and it is an upper Lipschitz continuous function. 

(H2) The curve S 2 is given and it is a lower Lipschitz continuous function. 

(H3) si(0) = a, 52 ( 0 ) = b, where a < b, and condition (iv) of problem ([T|) is not considered if a = 6. 
(H4) si(t) < S2{t)\/t € (0,T]. 

(H5) / is a non-negative continuous function in [a,b]. 

(H6) g and h are non-negative continuous functions in (0,T]. 


We will consider the following two regions: 

Dt = {{x,t)/si{t) <x < S 2 {t), 0 <t< T} and the so-called parabolic boundary = {(si(t),t),0 < 

t < T} U {(s 2 (t), i), 0 < t < T} U {(x, 0), a < a: < 6}. 

Definition 1. A function u is a solution of problem ([T]) ifu = u{x,t) verifies the conditions in ([T]) and 
1. u is defined in [ao,6o] x [0)7"], where uq = min{si(t), t G [0,r]} and b^ = max{s 2 (t), t G [0,T]}, 



2. u ^ CWdt = C{Dt) n W\{0,T)) n C^{Dt), where Vri((0,r)) = {/ G Ci((0,r]) : f G 

Li( 0 ,r)} 

3. u is continuous in Dt U dpDT except perhaps at (a, 0) and (6,0) where we will ask that 

0 < liminf u{x,t) < limsup u{x,t) < +oo 

(Xji)—>(a,0) (x,t)—^{a,0) 


0 < liminf u{x,t) < limsup u{x,t) < +oo. 

>-(6,0) (x,t)—>(6,0) 


Remark 1. We ask u to be defined in [ao,feo] x because the fractional derivative oDf‘u{x,t) 

involves the values of ufix^r) for all r in [0,t]. (See Figure 1). 


Remark 2. This kind of problems has not been studied in depth yet, but taking into account the results 
obtained in \2^ and \23\ . where some fractional Stefan problems has been solved explicitly, it is easy to 
check that the following problem 


oDfu{x, t) = Uxx{x, t) 0 < X < 0<f<T, 0<Q!<1 

tt(0, t) = B 0 < t <T 

u{t^/‘^,t) = C o<f<r 

admits the solution given by 


u{x, t) = B + 


C-B 




where Vk (•, — l) is the Wright function of parameters p = —% o-nd /3 = 1 defined by 


W{z;p-,fi) = J2 


k=0 


k\T{pk + fi)' 


z G C, p > -l,/3 G 


The function 1 — W{— ■, 1) is the ‘fractional error function”, so named because 


liml-fk(-x,-|,l)=er/(|) 


( 2 ) 


( 3 ) 


(see Theorem f.l). 


Hereinafter I will take A = 1, I will call to the fractional derivative in the Caputo’s sense of extreme 
0 = 0, oDf, and L" to the operator associated to the FDE 


:= 


d‘^ 

{dxY 




( 4 ) 


Proposition 1. If u is a function such that L°‘[u] > 0 in Dt, then u can not attain its maximum at 
Dt ■ 

Proof. Let us suppose that there exists (xo,to) G Dt (that is si(to) < 2:0 < S 2 {to), 0 < to < T), 
such that u attains its maximum at (a;o,to)- 

Due to the extremum principle for the Caputo derivative ( see |18)L we have that Dfu{xo,to) > 0. 

On the other hand, since u G C^{Dt), < 0. Then L“[ ^](^0) to) ^ 0 ) which is a contradiction. 

□ 


Corollary 1. If u is a function such that < 0 in Dt, then u can not attain its minimum in Dt- 








It is easy to adapt to the moving-boundary problem ([T]) results obtained in |16) for initial-boundary- 
value problems associated to the generalize FDE. For this reason we omit the proof of the following 
result. 

Theorem 1. Let u € CWd^ be a solution of O- Then either 

u{x,t) > OV {x,t) G Dt or u attains its negative minimum on dpDx- 
Let us enunciate the main result of this paper. 

Theorem 2. Let u G CWd^ be a solution of problem ([1]) satisfying the hypotheses {HI) — (HQ) 

1. If there exists to > 0 such that 

u{s 2 {to),to) = M = sup u, 

dpDT 

and, 

exists 6 > 0 such that |si(to) “ ' 82 (^ 0 )! ^ ^ u{x,to) < MVx G ( 52 (^ 0 ) “ b,S 2 {tQ)), 
then 

limint "("■*“) > 0. 

x/'s2{to) X — S2(to) 

If Ux exists at {s 2 {to),to), then 


( 5 ) 

( 6 ) 
( 7 ) 


Ux{s2{to),to) > 0 . 


( 8 ) 


2. If there exists to > 0 such that 


u{s 2 {to),to) = m = inf u, 
dpDr 


and. 


exists 5 > 0 such that |si(to) — •52(^0)! ^ ^ u{x,to) > mVx G (52(^0) “ b,S2{to)), 


then 


u{x,to) - u{s 2 {to),to) 

limsup- - —-- < 0. 

x/'s2{tQ) ® 'S2(^o) 


If Ux exists at {s 2 {to),to), then 

Ux{s2{to),to) < 0 . 

Proof. I will prove 1. The proof of 2 is analogous. 

Let us consider the function 

EaipiAt^) 


Wa{x, t) = e 


1 - exp{-/r(x - S 2 (to))}- 


+ M 


Haiti Af^) 

where A, p, and e will be determined and Ea is the Mittag-Leffler function defined by 

.,k 


Ea{z) = Y, 


k=0 


r{ak + l)' 


(9) 

( 10 ) 

( 11 ) 


( 12 ) 


(13) 


z G C, a > 0. 


(14) 








Note that Wa = M over the curve 


exp{-n{x - S 2 (to))} 


Ea {fJ-AE 


Observe that the curve (fT^ is the graphic of the function 


= 1 . 


f{t) = -ln 




+ 52(^0)) i S ( 0 ,to] 


^ V-®a(M^o). 

Clearly /(to) = S 2 {to) and it is easy to check that 

/ is an increasing function if /r > 0. 

Our next goal is to prove that there exists ti < to such that /(t) < S 2 (t) Vt G (ti,to). 


(15) 


(16) 


(17) 


Due to {H2), S 2 is a lower Lipschitz continuous function, then there exists a constant L > 0 such 


that 


Therefore 


S2(t) - S2(to) 
t - to 


<L V 0 < t < to- 


S2(t) > L{t - to) + S2{to) V 0 < t < to- 


(18) 


Taking into account that £'a(/rAt“) = X^^o r^ak+i) uniform convergent series over compact 

sets contained in (0,to] and that zT{z) = T{z + 1) Vz G C, we have that 

^ T[ak + 1) ^ T[ak + a) 

where the function E^^a is the generalized Mittag-Lefher function of parameters p = (3 = a dehned by 


EpA^) = 


k=0 


T{pk + (3) 


, z eC,p> 0,p eC. 


(19) 


Then, 


fit) 


1 1 pA 

fEaipAt^) tl-“ 


Ea,aAAn 


A Eg^aipAf^) 

t^-" EaipAf^) 


Now, let us define the function H : —?> M/ H{t) = 


( 20 ) 


77 is a positive function and it is a quotient of continuous functions, where its denominator is grater 
than 1 in [0,oo), then H is continuous in [0,oo). 

77(0) = > 0 because 0 < a < 1. 

77(+oo) = C > 0 because it is a quotient of continuous functions with equal order in 00 (see |10|i. 
Then, we can assure that there exists tuq > 0 such that 


H{t)>mo yt>0,\/A,p>0. (21) 

From p0]i and (1^ . we have that 

fito) > :^rno. 

h 

Selcting ^ > 0 such that ^4^7710 > L we can assure that 


( 22 ) 














fih) > L. (23) 

Lately, let be p > 0 / /'{to) — p > L. Due to the differentiability of / at to we can assure that there 
exists ti < to such that Vt G (ti,to), 

L < /'(to) - p < => /(t) < L{t - to) + /(to) = L{t - to) + S 2 (to) < S 2 (i) 

t — to 

Note that due to ([6|) and (fT7|l . we can select ti so that si(t) < /(t) < S 2 (t) Vt G (ti,to). 

Now, let be A(xi,to) (where xi = /(ti)), i?(s2(io))^o) and (^(xijti). Hypothesis (l6|) allows us to 
set ti again such that xi G ( 52 (^ 0 ) ~ ^ 2 (^ 0 )) and u < M m AC. 

Let be TZ the region limited by AB, AC and the portion of graph of / from B to C, which we will call 
CB.{See Figure 2) 

The region IZ^ = TZ° U {AB — {A, B}) and its parabolic boundary dpTZ = AC U CB will be con¬ 
sidered. 


Next, we will analyze the behavior of u and Wa in the parabolic boundary dpTZ. 

Let be Mq = max u{xi,t)- Because of the continuity of u, the hypothesis ([6|) and resetting ti if it is 

ti 


necessary, we can affirm that Mq < M. 

Calling p = 

M — Mq, it yields that 



u < M — 7] 

in 'AC. 

(24) 


u < M 

in CB. 

(25) 

By the other side. 

Wa = M 

in 

(26) 


In AC, considering that Ea[pAC) is an increasing function, we have that 

EaipAC) 


Wa{xi,t) = € 


1 - exp {-p(xi - S 2 (to))} 


EaipAt^) 


-I- M > e [1 - exp {-p(xi - S2(to))}] + M. (27) 


Taking e = ^ it results that 


From m we have that 


Wa{xi,t) > -r] + M. 


(28) 


D^{Ea{pAC)) = pAEaifiAC). (29) 

Now, applying the operator to the function Wa and using (|29j) it yields that 

L'^[wa]{x,t) = eexp{-p{x - S 2 ito))}^^p^p^{pA-<0 Ap = A + l. (30) 

Finally, we define the function z = Wa — urn. TZ. Let us analyze the behavior of z in the parabolic 
boundary dpTZ. 

From and (l28)l . 2 > 0 in AC. 

From (|25]) and (f26l) . 2 : > 0 in CB . 

Also, = L“[u)q] — L“[w] < 0 in TZi^. 








Applying Corollary [T| we can state that z cannot assume its minimum at ■ Then 

z >0 in TZ 

In particular, 

z{x,to) = Waix,to) - u{x,to) > 0, VXi < X < S2(^o)• (31) 

Recalling that u{s 2 {to),tQ) = Wais 2 {io,io)) = M, the next expression es equivalent to (I3ip : 


u{x,to) - u{s2{to),to) ^ Wa{x,to) - Wg(s2(tp, tp)) 
X - S 2 {to) ~ X- S 2 {to) 


Then 

u{x,to) - u{s2{to),to) ^ Wa{x,to) - Wa{s2ito, tp)) 

x/'s 2 (to) X — S 2 (to) x/'s 2 ito) X — S 2 {to) 

But Wa is a differentiable function at (s(to))io)) then 


Wa{x,to) - Wa{s2{to,to)) 
x/'s2{to) X — S2{to) 


{wa)xis 2 {to),to) = efj. = e{A + l) > 0 


and ([7|) holds. 


(32) 


Finally, if Ux exists at ( 52 (^ 0 ); ^o)) l|32p implies that 


Ux{s 2 {to),to) > {Wa)x{s 2 {to),to) > 0 (33) 

and dl]) holds. 

□ 

The same result is valid if we consider si instead of S 2 - 
Theorem 3. Let u € CWd^ ^6 ® solution of problem ([1]) satisfying the hypotheses {HI) — {H6). 


If there exists to > 0 such that 

u(si(to),to) = M = sup u, 

dpDT 

(34) 

and, 



exists (5 > 0 such that si(to) ■ 

- ^ 2 (^ 0 ) > hand u{x,to) < MMx G (si(to), si(to) + h), 

(35) 

then 

u{x,to) - u{si{to),to) 
limsup- - —-- < 0. 

(36) 

If Ux exists at (si(to)To); then 

Ux(si(to),to) < 0. 

(37) 

If there exists to > 0 such that 

u(si(to), to) = fTT-= inf u, 

dpDr 

(38) 


and, 


exists 6 > 0 such that \si{to) — S 2 {to)\ > 6 and u{x,to) > m\/x G {si{tQ), si{to) + 6), (39) 

then 

lirninf >0. (40) 

x/s-iito) X — Si(to) 

If Ux exists at (si(to),to)) then 


Ux{si{to),to) > 0 . 


( 41 ) 









3 An Application to Fractional Pree-Boundary Stefan Problems. 


Let us consider now the following fractional free-boundary Stefan problem for the FDE, where we have 
replaced the Stefan condition ^ = kux{s{t),t), t > 0, by the fractional Stefan condition 




's(t) = 

-kux{i 

5(t), t), t > 0, 0 < a < 1. 

(*) 

D°‘u{x,t) 

— '^XX 

(x,t) 

0 < X < s(t), 0<t<T, 0 <q;< 

{ii) 

u{x, 0) = 

fix) 


0 < X < b = s(0) 

{in) 

II 

o' 

9it) 


0 < t <T 

(iv) 

u(s(t),t) 

= 0 


0 <t <T 

[v] 

D“s(t) = 

kuxi 

s(t),t) 

0 < t < T, k > 0 constant 


(42) 


Definition 2. A pair {u, s} is a solution of problem (1421) if 

1. u and s satisfies (I42p . 

2. s is a continuous function in [0,T] such that s € VL(0, T) = {/ G (^^((O, T]) f f' ^ L^(0,T)}. 

3. u is defined in [0,6o] X [0,T] where 6o = max{s(t),0 <t< T}, 

U G CW£)rp . 

5. u is continuous in Dt U dpDx except perhaps at (0,0) and (6,0) where we will ask that 

0 < liminf u{x,t) < limsup u{x,t) < -|-oo 
>-( 0 , 0 ) >^( 0 , 0 ) 


and 


0 < liminf u{x,t) < limsup u{x,t) < -|-oo. 


6. There exists Ux{s{t),t) for all t G (0, T]. 

This kind of problems have been recently treated in (Sj IH [151 122])) and our goal now is to prove 
the next theorem involving the monotonicity of the free boundary. 


Theorem 4 . Let {ui,si} and {u 2 ,S 2 } be solutions of the fractional free-boundary Stef an problems 
corresponding, respectively, to the data {bi, fi, gi} and { 1 ^ 2 )/ 2 ) 52 }- Suppose that bi < b 2 , 0 < fi < f 2 
and 0 < gi <92- Then si{t) < S 2 (t)Vt G [0, T). 

Proof. Si and S2 are continuous functions and si( 0 ) = 61 < 62 = ■S2(0). Suppose that the set 
A = {t (z [ 0 ,T]/(si — S2){t) = 0 } 7^ 0 , and let be to = rnin^. Due to the continuity of si and S2, 
si{to) = S2(to)) and to is the first t for which si(to) = ■S2(to)- 

Let be h : [0, to] thn function h{t) = (si — S 2 ){t). h has the following properties: 

(h-l) h G C^(0, to] nClOjto] (due to definition [2|) . 

(h-2) h{0) = 61 — 62 < 0 . 

{h-3) /i is a non positive function and h{to) = 0. 

From (/i-l)-(/i-3), h attains its maximum value at to- 

Using the estimate ( eq. (12) of Theorem 2.1 of [2]), it results that 


Then, 


^ ^ “ t;^r(i - a) 

(43) 

D^h(to) > > 0. 

^ ^ “ t[^r(l - a) 

(44) 






Taking into account the linearity of the Caputo’s fractional derivative and that si and S 2 satisfies 
the Stefan condition (|42p -fvl. (|44p implies that 


U2xis2ito),to) - ui^{si{to),to) > 0 . 


(45) 


By the other hand, the function w{x,t) = U 2 {x,t) — ui{x,t) is a solution of the following moving¬ 
boundary problem 


D°^w{x,t) = Wxx(x,t) 
w{0,t) = {g 2 - gi){t) > 0 
w{si{t),t) = U 2 {si{t),t) 

_ w{x,0) = (/2 - fi){x) > 0 


0 < a; < si{t), 0 < t < to, 0 < a < 1, 
0 < t < to 
0 <t < to 
0 < X < bi = si(0) 


(46) 


Applying Theorem [T]to U 2 in the region where = {(x,t) / 0<t<to,0<x< S 2 {t)}, it results 
that U 2 {si{t),t) > 0. 

Now, applying again Theorem [T] to w in problem (1461) . it results that rc > 0 in where = 
{{x,t) / 0<t<to, 0<x< si(t)}. 

Then we can state that w attains a minimum at {si{to),to)- 

Now, if there exists e > 0 such that w{x,to) > 0 for all x € (si(to) “ £)Si(to))) applying Theorem [5}'2 
we can conclude that rc 3 ;(si(to); 0) < 0. And then 


U2x{s2{to),to) - uix{si{to),to) < 0, (47) 

which contradicts (j46p . 

If, by contrast, we have sequence {e^} such that e„ ^ 0 and for every re G N there exists Xn G 
(si(^o) “ en,si(io)) such that w{xn,to) = 0, then 

w{xn,to) - w{si{to),to) ^ ^ 

n-j-cxD 

and due to the existence of Wx{si{to),to) (definition [2]-6), it results that Wx{si{to),to) = 0. Then 


U 2 x{s 2 {to),to) - uix{si{to),to) = 0, (48) 

which contradicts (I46p again. 

This contradiction comes from assuming that there exists to > 0 such that to is the first t for which 
■si(to) = S 2 {to)- Therefore si{t) < S 2 (t)Vt G [0,T). □ 


4 Conclusions 

A one-dimensional moving-boundary problem for the time-fractional diffusion equation was presented, 
where the time-fractional derivative of order a G (0,1) was taken in the Caputo’s sense. Then, a 
generalization of the Hopf’s lemma, involving the behavior of the partial x-derivative of the solution 
at a boundary point, was proved. This last result was used to prove a monotonicity property of the 
free-boundary, when a free-boundary Stefan problem was considered. 
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